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Abstract 



According to the method, suggested in our previous work (arxiv:nlin. CD/0509012 ) and based 
on the consideration of the specially coupled systems, the possibility of physical realization of 
the phenomena of complex analytic dynamics (such as Mandelbrot and Julia sets) is discussed. 
It is shown, that unlike the case of discrete maps or differential systems with periodic driving, 
investigated in mentioned work, there are some difficulties in attempts to obtain the Mandelbrot 
set for the coupled autonomous continuous systems. A system of coupled autonomous Rdssler 
oscillators is considered as an example. 

1 Introduction 

It is known p!J |2], that complex analytic dynamics (CAD), studying behavior of complex maps, 
includes a lot of interesting phenomena, for example, presence of fractal Mandelbrot and Julia sets 
in the parameter and phase spaces. 

Let us start with a quadratic logistic map 

z'^X-z 2 , (1) 

where A is a complex parameter, and z is a complex variable. By definition, Mandelbrot set (fig. 1) 
is a set of points on a plane of complex parameter A, for which the orbit of an extremum z = 
of the map (JIJ during iteration procedure does not escape to infinity. The Mandelbrot set contains 
the so-called "Mandelbrot cactus" (designed on figure by gray color); this is a set of points in the 
parameter plane, for which the trajectory starting of the extremum of the map converges to a periodic 
attractor. 

"Mandelbrot cactus" consists of a big cardioid, corresponding to the existence of an attracting 
fixed point, and an infinite number of "leaves", corresponding to existence of attracting cycles of 
different periods. For example, the leaves of the doubled periods are placed along a real axis. 
The sequences of other period m-tupling (period-multiplication) bifurcations also can be found. In 
particular, the accumulation points for period-tripling and period-quadrupling bifurcation cascades 
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Figure 1: Mandelbrot set (a) and Julia sets for the quadratic complex map with different values 
of parameter: A = 0.5 (b), A = 0.8 (c), A = 1.42 (d), A = 0.5 + 0.7i (e), A = 0.123 - 0.745i (f), 
A = 0.0315 -0.7908i (g), A = -0.282 + 0.530i (h), A = 1.16 + 0.25i (i). The gray color designates the 
regions corresponding to existence of periodic dynamics (periods are marked by respective numbers); 
the black color designates points, at which the restricted in a phase space chaotic dynamics is 
implemented; the white color means the escaping of trajectories to infinity. 

have been investigated in the work of Golberg, Sinai and Khanin j^j. In the works of Cvitanovic and 
Mirheim |3J E] the universal properties of many other bifurcation cascades were studied. 

A bifurcation, which is responsible for originating a "leave", corresponds to a stability loss of a 
"parent" cycle, characterized by a complex multiplier with unit modulus and rational argument (in 
relation to 2tt). If the argument of a multiplier at the stability loss is irrational, then the domains with 
fractal boundaries, filled by invariant curves arise in the phase plane, the so-called Siegel disks [01 IB] • 

Fractal pattern close to the "Mandelbrot cactus" and denoted by black color in Fig. 1, corresponds 
to existence of the chaotic dynamical regimes in the phase space. 

In figures l(b-i) the Julia sets for different values of complex parameter A are shown. The Julia 
set is a border between basins of attraction to infinity (white color) and to a periodic motion (gray 
color) in a plane of complex variable z. One can distinguish the following types of Julia sets: 

• for values A, belonging to the "Mandelbrot cactus", the Julia set is connected and enclose an 
interior basin (figs. l(b,c,f-i); 

• for values A, at which chaotic dynamics exists, the Julia set is also connected, but has no inner 
region (fig. Id); 

• for values A, outside the Mandelbrot set, the Julia set is disconnected (fig. 1/). 

It is obvious, that ID complex map can be represented equivalently by a 2D real map (for this 
purpose it is necessary only to separate real and imaginary parts of the equation). However, the 
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mentioned phenomena of CAD are intrinsic only to a very special class of the real 2D maps, namely 
for the analytic maps, obeying the Cauchy-Riemann conditions. Violation of the analyticity leads to 
drastic changes of the dynamics of the map [OJ HOI CTJ [T21 HB! - Thereby, a following problem arises: Is it 
possible to specify actual physical systems demonstrating phenomena of CAD? Recently, this problem 
attracts great attention. The physical applications of complex dynamics for such problems, as the 
renormalization group approach in the theory of phase transitions and the theory of a percolation 
were discussed [HI HH1 EH Ej • In the paper of Beck [THj a theoretical possibility of the construction 
of the physical system, in which the Mandelbrot set would arise, was considered. The suggested 
approach is based on analysis of a motion of a charged particle in a double-peak potential with 
non-linear damping. The particle is driven by magnetic field, depending on time and on the particle 
velocity, and effected by external shot pulses, time intervals between which also depend on the particle 
velocity. 

In works [T5J 120] and in the present work, we offer a simpler and universal approach of constructing 
models manifesting the Mandelbrot set and other phenomena of CAD, which may be designed as 
realistic physical systems. This method allow us to carry out a physical experiment and present the 
first observation of the Mandelbrot set (21] . The special structure of the Fourier spectrum of signal 
generated by experimental system at the period-tripling accumulation point is presented in [2"2"j . 

The method developed in [TH] and in the present work is based on using coupled systems, demon- 
strating transition to chaos through period-doublings. It is known that such behavior is characteristic 
for a very wide class of nonlinear dissipative systems of various physical nature. 

As shown in ^Hj, it easy enough to arrange realization of the Mandelbrot set for systems with 
discrete time and non-autonomous periodically-driven systems (coupled logistic and Henon maps and 
coupled oscillators with quadratic nonlinearity and harmonic external driving have been investigated). 
For this one considers two identical elements demonstrating the period-doublings with symmetrical 
coupling, arising from complexification of the variables and the control parameter, responsible for the 
period doubling in the original system. Domain of generalized partial synchronization (defined by 
the term, that the dynamical variables of subsystems do not escape far from each other) corresponds 
to the Mandelbrot set. 

Such special kind of coupling provides a special symmetry in the multi-dimensional system, which 
is necessary for implementation of the analyticity conditions for the discrete time model map (or the 
stroboscopic Poincare map for the continuous system). It is a simple problem to construct the system 
with such coupling in comparison to the system, suggested by Beck. 

From our point of view, the problem of realization of phenomena of CAD for autonomous con- 
tinuous systems is now opened and seems interesting. Consideration of this problem is a main goal 
of the present paper. In section 2 we reproduce the method of the obtaining of the coupling function 
for the logistic maps. In section 3 we try to implement a similar procedure to a continuous flow 
system, namely, to the Rossler oscillators. 



3 



2 From complex quadratic map to coupled logistic maps 



Let us start with the notion that one-dimensional complex quadratic map is equivalent to the system 
of two real coupled quadratic maps with a special type of coupling. 

Separation of real and imaginary parts in the complex equation JIJ yields 

Z re > \ re Z re -\- Z^ m: ^im * Aj m ^^re^im- (2) 

Next, we introduce the following designations 

X\ z re -\- (3zi m} %i z re f3zi m} 

Al = A re + f3\mi ^2 = A re — /3Aj m . 

As a result we obtain a system of two coupled quadratic maps 

x'i -»• Ai -x\ + e(x 2 -xi) 2 , ^ 
x' 2 ^ \ 2 - x\ + e(xi - X2) 2 , 

where e = (1 + /3 2 ) / 4/3 2 is the parameter of coupling. Note a special type of coupling in these 
equations. It can be interpreted as an identical simultaneous shift of control parameters in both 
partial systems, proportional to the squared difference of dynamic variables at each step of discrete 
time. 

It is worth nothing that the coefficient e = (1 + /5 2 ) /4/3 2 for any f3 is larger than 1/4. Nevertheless, 
formally we can investigate system (jlj) with any e. 

In fig. 2 we present the charts of the parameter plane (Ai, A2) for the coupled maps (jlj at several 
values of parameter e. One can see the usual Mandelbrot set, rotated by 45°, takes place at e = 0.5. 
For 0.25 < e < +00 we have a distorted Mandelbrot set on the parameter plane. The cactus leaves of 
this Mandelbrot set correspond to existence of periodic motion of different periods. At e — 0.25, the 
set on the parameter plane, for which the point starting from the origin does not escape to infinity, 
looks like a set of strips, where the period doubling cycles occur. At e < 0.25 it transforms to a 
rhombus-like structure. At a particular e = (uncoupled logistic maps) it is a square. 

The generalization of coupling to the case of —00 < e < +00 corresponds to the original map ((TJ), 
variable and parameter of which are so-called two-component numbers [231 El 123 1213 123 EH] • This 
is a special algebraic system, which elements are defined as follows 

z = x + iy, i 2 = a + ib, where a, b G R. (5) 

According to j2Hl> there are three special cases: i 2 = — 1 - the usual complex numbers, i 2 = +1 
- the so-called perplex numbers, i 2 = - the dual numbers. All other algebraic number systems are 
isomorphic to complex, perplex or dual numbers depending on, whether the value of (a + b 2 )/4b 2 
is positive, negative or zero, and are known as elliptic, hyperbolic or parabolic number system, 
respectively. In terms of parameter e these conditions look as follows: 

1) the case e > 0.25 corresponds to elliptic numbers isomorphic to complex numbers, implemented 
at e — 0.5; 
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Figure 2: The charts of the parameter plane (Ai, A2) for the coupled logistic maps (J3J) with different 
values of parameter of coupling: e = —0.1 (a), e = 0.0 (b), e = 0.1 (c), e = 0.25 (d), e = 0.3 (e), 
e = 0.5 (f), e — 1.0 (g). The figures (a-c) correspond to hyperbolic numbers, figure (d) - to parabolic, 
and figures (e-g) - to elliptic numbers. 
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2) the case e < 0.25 corresponds to hyperbolic numbers isomorphic to perplex numbers, implemented 

at e = 0; 

3) the case e = 0.25 corresponds to parabolic or dual numbers. 

Thus, the existence of three topologically different structures on the plane of parameters (Ai, A2), 
namely, fractal structure, similar to Mandelbrot set, rhombus-like structure and system of strips, is 
explained by existence of three different algebraic systems of numbers. 



3 Coupled autonomous Rossler oscillators 

Let us consider autonomous flow system — Rossler oscillator 



x = -(y + z), 

y = x + ay, (6) 
z = b + z(x — c), 



where x, y, z - are dynamical variables, a, b, c - are parameters [29]. It is well-known fact, that 
for the map at the Poincare cross-section for this system, defined for example by the plane y = 0, 
the transition to chaos through period-doubling bifurcations is possible. At figure |3] the charts of 
the parameter (c, a) (a) and (c, b) (b) planes are shown. Cascade of period-doubling bifurcations 
is visible, for example, with changing of parameter c with fixed a and b. Let us take for example 
the values of parameters a = 0.2, b = 0.2. The bifurcation tree for this values of parameters are 
represented at Fig. El (c). 

Let us construct the system of coupled Rossler oscillators. To obtain the function of coupling we 
use the procedure described in previous section. The complexification of Rossler system (in such way, 
that dynamical variables x, y, z and parameter, responsible for the period doublings c are complex, 
and other parameters a and b are real) with introduction of the variables and parameters designations 
like (JHJ) gives the coupled systems of following form 

xi = -(yi + zi), x 2 = -(y 2 + z 2 ), 

iji = x x + ay 1 , y 2 = x 2 + ay 2 , 

zi=b + zi(x\ - ci)— z 2 = b + z 2 (x 2 - c 2 ) — 

-e{z 2 - zi)((x 2 - xx) - (02 - ci)), -e(zx - z 2 )((xx - x 2 ) - {c x - c 2 )). 

But the attempt to obtain the Mandelbrot set at the plane of new parameters (cx,c 2 ) of the 
system (JZj) has failed. The investigation of parameter space has not given suitable results - the 
domain of generalized synchronization has appeared not to be similar to Mandelbrot set and is 
represented by thin band alone the diagonal line c x = c 2 (this band become more and more thin 
with coupling parameter growth). Apparently, the reason of such behavior lie in the existence of a 
detuning of the phases of coupled subsystems. In following subsections we try to avoid the problem 
of phase detuning by several special approaches for construction of the coupled systems. 
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Figure 3: Charts of the parameter plane (c, a) with b = 0.2 (a) and (c, b) with a = 0.2 (b) for the 
Rossler oscillator. At the fig. (c) the bifurcation tree of the Rossler oscillator at the plane (X n , c) is 
represented. The value X n is the discrete stroboscopic cross-section, corresponded to the values of 
dynamic variable x at the Poincare section y = with x < 0. 



3.1 Truncated system of coupled Rossler oscillators 

To investigate the synchronization regimes of the coupled Rossler systems often it is more convenient 
to use a cylindrical system of coordinates [3U], in which the set of variables (x,y,z) exchanges by 
variables A, cp, z, where <p = arctan ^ - is a phase, A = (x 2 + y 2 ) 1 ^ 2 - is an amplitude. The system © 
in such cylindrical coordinates is represented as follows: 

A = aA sin 2 <p — z cos <p, 

ip = 1 + sin (p cos <p + 4 sin (p, (8) 
z = b — cz + Az cos <p. 

Poincare cross-section in this case is defined as <p = 2im, n = 1, 2, ... 

The coupled Rossler oscillators (jZj) rewritten in the cylindrical coordinates looks as 



A\ = aAi sin 2 <p% — z\ cos <p\, A 2 = aA 2 sin 2 <p 2 — z 2 cos (p 2 , 

Zi = b + z 1 (A 1 cos^i - ci)- z 2 = b + z 2 (A 2 cos<p 2 - c 2 )- 

-e(z 2 - z 1 )((A 2 - Ax) cos^i - (c 2 - ci)), -e{z\ - z 2 ){{A x - A 2 ) cosp 2 - (c x - c 2 )), 

(p\ = 1 + a sin <px cos (pi + sin (p\ (p 2 = 1 + a sin <p 2 cos (p 2 + sin (p 2 . 



(9) 



Let us simplify this system, neglecting the second and third terms at the equations for phases (Such 
approach does not connect with some kind of physical interpretation and is considered exclusively 
for the showing up the relation between phase synchronization and possibility of Mandelbrot set 
realization). In this case, we have only one variable of phase depending on time by the linear law 
<p = 1. Then the initial system converts to the form, similar to nonlinear system with external 
harmonic driving. Thus, the simplified system looks as 

A\ = aA\ sin 2 (p — z\ cos <p, A 2 = aA 2 sin 2 (p — z 2 cos <p, 

z\ — b + zi(Aicos(^ - ci)— z 2 = b + z 2 (A 2 cos ip - c 2 )— 

-e(z 2 - z t )({A 2 - Ax) cos <p - (c 2 - c x )), -e{z x - z 2 )({A 1 - A 2 ) cos<p - (c x - c 2 )), 

= 1. 
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Figure 4: Chart of the parameter plane (ci, c 2 ) for truncated system of coupled Rossler oscillators 
(fTUj) with e = 0.5, a = 0.2, 6 = 0.2. 




X", 






(a) 



(b) 



(d) 



Figure 5: Basins of attraction at the phase plane (x\ = A\ cos<^, x 2 = A 2 cosy?) for truncated system 
of coupled Rossler oscillators (jlO)) with e = 0.5, a = 0.2, b = 0.2 for the values of parameters c\ 
and c 2 corresponding to existence of an attracting fixed point at c\ = c 2 = 1.0 (a), cycle of period 
2 at C\ = c 2 = 2.25 (b), cycle of period 4 at c\ = c 2 = 2.58 (c) and cycle of period 3 at c\ = 0.75, 
c 2 = 2.25 (d) in Poincare cross section. 
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(d) 

Figure 6: Phase portraits for the coupled Rossler oscillators (|llj) with the same values of parameters 
as for fig. El The attractors of period 1 (a), 2 (b), 3 (c,d) and 2 (e, f) are represented. 
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Let us term it as truncated system of coupled Rossler oscillators. It is known from the work ^H] , that 
for such systems the phenomena of CAD can be realized. Really, the numerical simulations of (fTU|) 
show, that at a plane of parameters (01,02) (fig. 4) one can see the structure similar to Mandelbrot 
set, and at the plane (x±, x 2 ) (fig- 5) - structures, similar to Julia sets. The present fact confirms the 
guess of necessity of the phase synchronization of coupled subsystems. At figure 6 the projections of 
several periodic attractors (including period-three attractor) in the phase space are shown. 

3.2 Unidirectionally coupled Rossler systems 

Now, let us consider the system of coupled Rossler oscillators with the mutual equation for the phase 

A\ = aA\ sin 2 y? — zi cos ip, A 2 = aA 2 sin 2 y? — z 2 cos y?, 

zi = b + z 1 (A 1 cosy? - c x )- z 2 = b + z 2 (A 2 cosy? - c 2 )- 

-£(2:2 - z 1 )((A 2 - Ax) cosy? - (c 2 - ci)), -e(zi - 2 2 )((^i - A 2 ) cosy? - (ci - c 2 )), 

y? = 1 + sin ip cos y? + j^- sin ip. 

The similar systems consisting of two coupled subsystems and the equations, mutual for these subsys- 
tems were explored in the works of Pecora and Caroll jHU E21 ESj • In these works the N-dimensional 
dynamical system u = f(u) is subdivided on two subsystems {u} = {v, w}, one of which is duplicated 
as follows v = g(v,w), w = h(v,w), w' = h(v,w'). Then, it is possible to view {w'} as response 
slave system drove by a signal from master system {v,w}. The major interest represents the study 
of synchronization regimes of this system with a driving signal jSUESj- 

The numerical experiment show, that the domain of generalized synchronization for system (fT2|) 
has some features, appropriate for CAD. For example, one can see the leave of period 3 at the 
parameter plane (fig. 6). The leaves of higher periods are small enough and are not visible on the 
represented chart. 

At fig. 7 the basins of attraction at a cross-section plane of a phase space of system x\ = A\ cosy?, 
x 2 = A 2 cos y? are represented. Non-symmetrical structure of the domain of generalized synchro- 
nization at fig. 6 and the basins of attraction at fig. 7, can be explained by unidirectional nature of 
coupling between driving and response subsystems. 

From a point of view of the experimental physical applications, dynamics of the system (fT2|) 
requires further study. 
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Figure 7: Chart of the parameter plane (c\,C2) for coupled Rossler oscillators (fTTj) with e = 0.5, 
a = 0.2, b = 0.2. 
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Figure 8: Basins of attraction at the phase plane (x\ = A\ cos<£>, %2 = A2 cos(p) for the system (fTTJ) 
with e = 0.5, a = 0.2, b = 0.2. The values of parameters c\ and c<i corresponds to existence of an 
attracting fixed point with c\ = ci = 1.5 (a), cycle of period 2 with c\ = C2 = 3.42 (b), cycle of 
period 3 with c\ = 2.92, c 2 = 0.83 (c) and cycle of period 2 at symmetric point C\ = 0.83, c 2 = 2.92 
(d). 
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Figure 9: Phase portraits for the coupled Rossler oscillators (fTTj) with the same values of parameters 
as for fig. El The attractors of period 1 (a), 2 (b), 3 (c,d) and 2 (e, f) are represented. 
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3.3 Coupled Rossler oscillators with additional coupling, which synchro- 
nize phases 

The next our step is to consider two different equations for phases (fx and ip 2 for each of the partial 
systems. To solve the problem of phase detuning we enter the additional coupling with parameter 
5 to the equations for phases, which thus, provide phase synchronization of Rossler subsystems and 
does not influence amplitudes 

A\ = aA 1 sin 2 ip 1 — z\ cos ipx, A 2 = aA 2 sin 2 (p 2 — z 2 cos (p 2 , 

z 1 — b+ z\(Ai cost^i - ci)— z 2 = b + z 2 (A 2 cos^ 2 - c 2 )— 

-e(z 2 - zt)((A 2 - Ax) cosy?i - (c 2 - c x )), -e{z\ - z 2 )((Ax - A 2 ) cosip 2 - (ci - c 2 )), (12) 

(px = 1 + a sin <£>i cos ^ + sin ^ 2 = 1 + a sin yj 2 cos <^ 2 + ^ sin <^ 2 

+5sin(v? 2 - (fi), +5sin(ipi - <p 2 )- 

In usual rectangular coordinate system {x = Acosip, y = A sirup, z) the equations (fT^j) looks as 
following 

x x = -(yx + zx) - 5yx {x s + ^2 { lf_l y 2 )lf i, x 2 = -(y 2 + z 2 ) - 5y 2 ^^r^^m > 
y x = xi + ayx + Sxx j^^ff^ , h = x 2 + ay 2 + 6x 2 {x 2 + ^}^ y 2 )1/2 , ^ 
zx = b + zx(xx - cx)— z 2 = b + z 2 (x 2 - c 2 )- 

e{z 2 - zx){{x 2 - xx) - (c 2 - cx)), e(zx - z 2 )((xx - x x ) - (c x - c 2 )). 

The numerical simulations show, that at a plane (ci, c 2 ) for system (|T^|) the arising of fractal set, 
similar to Mandelbrot set (fig. 8) is possible. Nevertheless, it is easy to see the considerable distortion 
of its configuration. For explanation of this fact we may assume, that the entering of additional 
coupling breaks the conditions of complex analyticity for the resulting stroboscopic Poincare map. 

At fig. 9 the basins of attractions similar to Julia sets are represented. At fig. 10 the attractors of 
period 1 (a), 2 (b), 4 (c) and 6 (d, e) are shown. Let us note, that with identical values of parameters 
Cx and c 2 (see fig. a-c) the values of variables in partial systems are coincided, that corresponds 
to the full synchronization. In a case ci 7^ c 2 (fig. d, e), when the phenomena of CAD can be 
realized, the values of dynamical variables of partial systems do not coincide, that corresponds to 
realization of generalized partial synchronization. Thus, the phenomena of CAD such as period- 
tripling bifurcations can be implemented, when the point (ci, c 2 ) belongs to the region of generalized 
partial synchronization. 
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Figure 10: Chart of the parameter plane (ci, C2) and its enlarged fragment demonstrating self- 
similarity for coupled Rossler oscillators (fT3j) with e = 0.5, a = 0.2, 6 = 0.2, 5 = 1.0. 
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Figure 11: Basins of attraction at the phase plane (xi,x 2 ) for the coupled Rossler oscillators (|T3|) 
with e = 0.5, a = 0.2, 6 = 0. 2, 5 = 1.0 for different values of parameters c\ and C2 corresponding to 
existence of an attracting fixed point at C\ = C2 = 2.0 (a), cycle of period 2 at c\ = c<i = 3.5 (b), 
cycle of period 4 at c\ = C2 = 4.0 (c) and cycle of period 6 at c\ = 2.98, C2 = 3.9 (d). 




(d) (e) 

Figure 12: Phase portraits for the coupled Rossler oscillators (fT3|) with the same values of parameters 
as for fig. 13. The attractors of period 1 (a), 2 (b), 4 (c) and 6 (d, e) are represented. 
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3.4 Coupled Rossler oscillators with different time variables 

Let us suggest one more method of obtaining of Mandelbrot set for autonomous system with continu- 
ous time, such as the coupled Rossler oscillators, which allow to avoid the problem of phase detuning 
and doesn't break analyticity conditions. This method lies in following: let us complexify Rossler 
system by another way. Let us consider the variables A, z and time t as complex and the variable ip 
as real numbers. For this purpose, it is suitable to rewrite the original Rossler system in cylindrical 
coordinates as 



dA aA sin i£— z cos ip_ 

dp 1+sin ip cos y+Jf sin ip ' 

dz b— cz+Az cos ip 

dp 1+sin ip cos tp+ sin ip ' 
dt 1 



(14) 



(15) 



dp 1+sin p cos y+ sin ip ' 

After complexification and the variables and parameter designations we obtain the system of coupled 
reorganized Rossler oscillators 

dAi _ Ai (aAi sin 2 y— z\ cos ip) . dA 2 _ A 2 (aA 2 sin 2 ip~z 2 cos ip) , 
dp Ai+A\ sin ycos y+zi sin y dp A 2 +A 2 sin p cos p+z 2 sin y 

+ef 1 (A 1 , A 2 , zi, z 2 , ci, c 2 ), +ef 1 (A 2 , A t , z 2 , z 1 , c 2 , ci), 

dzi _ Ai(b-zi(ci-Ai cosy)) . dz 2 _ A 2 (b~z 2 (c 2 +A 2 cos p)) . 
dp Ai+Ai sin p cos y+zi sin ip dp A 2 +A 2 sin p cos p+z 2 sin y 

+ef 2 (Ai, A 2 , zi, z 2 , ci, c 2 ), +ef 2 (A 2 , A u z 2 , z ± , c 2 , ci), 

dti Ai I dt 2 Ag I 

dip Ai+Ai sin ip cos ip+z\ siny dp A 2 +A 2 sin ip cos p+z 2 sin tp 

+ef 3 (A 1 , A 2 , z lt z 2 , ci, c 2 ), +ef 3 (A 2 , A u z 2 , z 1 , c 2 , ci). 
where fi, f 2 and f 3 are the functions of coupling, which can be expressed as follows 

/l(^l,2, A 2t i, Z lt2 , Z2,l, Cl, 2 , C 2j i) = 

~2u ly2 u 2yl (A lt2 (w li2 -w 2il )+A 2yl (w ly2 +w 2yl ))+2ul 2 A 2tl (w li2 -w 2il )+2u^ 1 A lt2 (w lt2 +w 2tl ) 

«i,2(«f i2 + u 2,i) ' 

f2{A lt2 ,A 2 

,1, Z\,2, Z2,l, Cl,2, C 2 ,lJ — 
^-2u ia u 2tl (A li2 +A 2}1 )+2ul 2 A 2A +2u^ 1 A h2 

u iM u i, 2 + u ' 2 ,i) 

I -"i,2M2,i(Ai i2 (-zi,2((ci,2-^i,2 cos <p)-(c 2tl -A 2 i cosy))+2 2 ,i((c li 2"Ai i2 cosy)+(c 2 ,i-A 2 ,i cosy)))) 

«l,2("i j2 + u 2,l) 

I -"i,2«2,i(A2,i(-zi,2((ci,2-^i,2 cosy)+(c 2 ,i-A 2 i cosy))+z 2 ,i((c li 2-A li 2 cosy)-(c 2 ,i-A 2 ,i cosy?)))) 

^ »y(«i 2 +«y ( 16 ) 

« 2 2 A 2 ,i (zi,2((ci,2 -Ai, 2 cos y)-(c 2 ,i -A 2 ,i cos y))-z 2 ,i ((ci, 2 -Ai j2 cos y)+(c 2 ,i -A 2 ,i cos p))) 

H ' "i,2K, 2 +^,i) h 

uf 2 Ai >2 (-2i j2 ((ci, 2 — A i>2 cosy)-(c2,i-A 2 ,i cos y))+z 2 ,i ((ci >2 — A i>2 cos y)+(c 2 ,i-A 2 ,i cosy))) 

Ml,2(^ i2 + U 2,l) 

2u\ 1 A\ >2 z\. 2 (c\. 2 -A\ t2cas v ) 

ui, 2 ( u l, 2 + u l,i) ' 

,1; ^1,2, ^2,1, Cl )2 , C 2; iJ — 
-2ui i 2M2,i(Ai,2+A2,i)+2« 2 i2 A 2 ,i+2u 2 . j1 Ai j 2 

«l,2(«f i2 + u 2,l) 

where 



(17) 



wi,2 = ^1,2 + ^4i,2 sin y? cos </? + ^i, 2 sin </?, 
u>i j2 = aA 12 sin 2 </? — ^ lj2 cos ip 

The represented unusual way of complexification appear to be successful by some reasons. At first, 
let us remind that we miss the problem of phase detuning. Thus, the variable ip responsible to the 
Poincare cross-section {ip = 2im, n = 1,2,...) remains real (it is not suitable to consider the variable 
y as real, i.e. to consider the derivatives dx/dy, dz/dy, dt/dy, because y is an non-monotonous 
function of time). Second, such complexification does not destroy analyticity. 
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Figure 13: The chart of the parameter plane (ci,c 2 ) for the coupled Rossler oscillators ([T5|). 
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Figure 14: Basins of attraction on the phase planes (xi,X2) for the coupled Rossler oscillators (fTKjl . 
The basins of attractors of period 1 with c\ = c 2 = 1.7 (a), 2 with ci = c 2 = 3.35 (b),4 with 
Ci = c 2 = 4.0 (c) and 3 with ci = 3.222, c 2 = 0.828 (d) are represented. 
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Figure 15: Phase portraits for the coupled Rossler oscillators (JT5j) with the same values of parameters 
as for fig. 16. The attractors of period 1 (a), 2 (b), 4 (c) and 3 (d,e) are represented. 
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Figure 16: The dependence of time detuning on versus the phase with the values of parameters 
Ci = 3.222, c 2 = 0.838 corresponded to existence of 3-cycle. 

The represented at fig. 11 chart of the parameter plane (ci,c 2 ) demonstrates the full agreement 
with usual Mandelbrot set. The basins of attraction of cycle of period 1 (fig. 12a), 2 (fig. 12b), 3 
(fig. 12c) and 9 (fig. 12d) are also similar with Julia sets. The next figure 13 shows the attractors for 
the subsystems represented in rectangular coordinates. Note, that if C\ ^ c 2 , the variables of time 
for each subsystem have no coincided values, i.e. the trajectories of A\,z\,tp and A 2 ,z 2 ,(p run up 
to the Poincare cross-section <p = 27m not simultaneously. Fig. 14 demonstrates the dependence of 
time detuning on the variation of phase with the values of parameters ci, c 2 corresponded to cycle of 
period 3. As it is visible, the value of t 2 —t\ oscillates periodically near the averaged value increasing 
by linear law. 

Also, it is necessary to note, that the variable of time are not subsisted in the equations for 
amplitude A and variable z, and doesn't influence its dynamics. Therefore, the equations for t\ and 
t 2 can be rejected. Thus, system (JT5)l is also the truncated system. Besides, unfortunately, it is 
difficult to realize the system (|15jl experimentally. 

4 Conclusion 

In present work we have made an attempt to obtain the phenomena of CAD such as fractal Mandel- 
brot and Julia sets in the autonomous continuous systems. As an example the autonomous Rossler 
oscillator has been considered. 

It is shown, that complexified system can be represented as two coupled real systems. Such 
representation is useful for some reasons. At first, it simplify construction of a real physical model. 
Secondly, with entering of special coupling to the system of two identical devices of any nature demon- 
strating the period-doubling cascade, one may expect the whole system to demonstrate phenomena 
of complex analytic dynamics. 

The main result of the work is that the realization of the Mandelbrot set in parameter space 
of the autonomous continuous systems requires provision of additional conditions such as phase 
synchronization. The phase synchronization can be obtained by the truncation of the system of 
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equations, by entering of phase coupling or by the special type of complexification. 

Nevertheless, the problem of the realization of Mandelbrot set in autonomous flow systems require 
more careful study. For example, one can assume, that the codimension of phenomena relevant to 
CAD is higher for the flow systems and Mandelbrot set as domain of generalized synchronization 
can be observed not on a plane, but on some surface in parameters space. It is necessary to study 
in more detail the regimes of synchronization for the offered system and the process of influence of 
a phase detuning on structure of parameter space. 
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